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ABSTRACT 

We present the results of a search for occultation events by objects at distances 
between 100 and 1000 AU in lightcurves from the Taiwanese-American Occultation 
Survey (TAOS). We searched for consecutive, shallow flux reductions in the stellar 
lightcurves obtained by our survey between 7 February 2005 and 31 December 2006 
with a total of ~ 4.5 x 10^ three-telescope simultaneous photometric measurements. 
No events were detected, allowing us to set upper limits on the number density as a 
function of size and distance of objects in Sedna-like orbits, using simple models. 
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Introduction 



During and just after the epoch of planet formation, a significant number of planetesimals 
were left that had not accreted into planets. Some of the planetesimals near the giant planets 
were ejected to orbits with large semi-major axes, and many more were hurled hyperbolically into 
interstellar space. Those objects which formed beyond p roto-Neptune would probably remain near 
the ecliptic plane a nd became part of the Kuiper Belt (jOorg ll950|; lKuipeiill974l : lFernandeall98Cl : 
Duncan et al]|l987l ). 

Up to the present, more than 1000 Kuiper Belt Objects (KBOs) have been discovereclll. The 
estimated total mass of the Kuiper Belt from observations to date is ~ 0.1 M^, which is a few orders 



of ni agnitude less than the mass predicted for the minimum mass solar nebula (jWeidenschilling 
19771 ). Theoretical simulations show that in order to assemble 1000- km size objects between 30 
and 5 AU, initia lly there must have been much more mass in this region ([Kenyon &: Bromley 
2004 : ISternI l2005l ) tha n the current observat i ons suggest. Many surveys indicate a Kuiper Belt 
"cliff" beyond 50 AU (jchiang fc Brownlll99d : lAllen et ahlboOll : bernstein et al.ll2004l ). Given that 
more than 99% of the original disk mass appears to have been depleted, some sort of dynamical 
perturbations must have removed objects from this region. 

Furthermore, t he origin of 2003 VB 12 (Sedna, with orbital parameters of a=531 AU, q=76 AU, 
i = 12°) is puzzling (iBrown et al.ll2004al ) . Searches for thermal radiation from Sedna with the IRAM 
30 m telescope and Spitzer have put an upper limit on its diameter of 1600 km (jBrown et al.ll2004bl ) . 
The mass of Sedna, from the upper limit on the size and mean density of KBOs, is estimated to 
be < 10~^ M3. Given its large size and highly eccentric orbit, in situ formation seems unlikely. 
Sedna is unique because its perihelion at 76 AU is far outside the reach of dynamical perturbation 
by Neptune, and the aphelion of ~1000 AU is also too small for Sedna to be affected sufficiently by 
giant molecula.r clouds or galactic tides in the solar neighborhood. Meanwhile, recent observations 
([Barucci et al.ll2005l : lEmerv et al.ll2007l ) have shown that Sedna is likely to have meth ane, water and 



nitro gen ices. The similarity of the spectroscopic features between Sedna and Triton ([Barucci et al 
2OO5I ) suggests its formation near the giant planet region. 



Several theories have been pro posed to explain the existence of Sedna in its present orbit. For 
instance, a close stellar encounter (jFernandez &: Bruninill2000l : llda et al.ll200d : iKenvon &: Bromlev 
2004 : iMorbidelli &: LevisonI 12004 ) could have depleted the Kuiper Belt and excited outer Solar 
System objects, such as Sedna, to their present orbits. Other possibilities have also been sug- 
gested, including rogue mars-sized bodies, or multiple stellar encounters in the Sun's birth c luster 
(jBrown et al]l2004al : iMorbidelh k LevisonI I2OO4I : loiadman fc Chanlbood : iBrasser et aDl2007l ). 



Note that all of these theories also predict that many more objects would be found in Sedna-like 
orbits, i.e. orbits with perihelia large enough to be unaffected by Neptune and aphelia small enough 



'^See lhttp: //www. cf a.hcLrvcLrd. edu/iau/lists/THOs .html| for a list of these objects. 
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to be unaffected by galactic tidal forces or giant molecular clouds. Furthermore, the fact that Sedna 
was discovered v e ry clos e to perihelion, where it spends a very small fraction of its orbital period, 
led iBrown et al.l (|2004al ) to infer that there may be as much as many as 500 objects in similar 
orbits with the same size as Sedn a or larger. However, recent search attempts ([Larsen et al.l 120071 : 
Brownll2008l: ISchwamb et al.ll2009l ) did not find any new member of this population. In their recent 



survey. 



Schwamb et al 



(120091 ) modeled objects with the same semi-major axis and eccentricity as 
Sedna. A best fit of 40 objects brighter than or equal to Sedna in that region is consistent with 
their results of detecting one object with a perihelion past 70 AU. 

Most objects in Sedna-like orbits would have escaped from direct detection due to their enor- 
mous distances from the Earth. Direct detection of this kind of object is limited to relatively 
close objects and by large telescopes, so the m ajority of the population is still well beyond the 
detection limit for most ground-based telescopes. iBaileyl (jl976l ) suggested that the existence of the 
"invisible" objects in this region might be inferred by stellar occultation, and the utility of using 
robotic telescopes to search for occulta tions by outer Solar System objects has been demonstrated 



(Axelrod et al 



19921 : lCooketalJll995|). 



The Taiwanese- Ameri can Occultation Su r vey (TAOS) aims t o investigate the size distribution 



of small (~1 km) KBOs (ILehner et al 



200S 



Zhang et al.l 120081 ). Four 50 cm robotic telescopes 



have been installed at Lulin mountain in central Taiwan. Each telescope is a fast F/2 system with 
a 3 deg'^ field of view, equipped with a 2kx2k CCD camera that can perform 5 Hz photometric 
sampling on a few hundred stars simultaneously. Multiple telescopes are used to reduce the false 
positive event rate. The TAOS system has been in routine operation with three telescopes since 
February 2005, and four-telescope observations began in August 2008. 

TAOS was designed to search for KBOs near 50 AU, for which a typical occultation by a KBO 
of a few km across would manifest itself as a flux reduction in only one or two consecutive time 
series measurements. With a nu ll detection, a strin gent constraint has been set on the number 
and size distribution of KBOs by IZhang et al.l (120081 ). However, the same data set can be used to 
detect objects in the distant Solar System all the way to ~1000 AU; i.e., the survey is also sensitive 
to objects in Sedna-like orbits. For objects at 100 - 1000 AU, the Fresnel scale and the projected 
size of the background star a re larger than for objects in the Kuiper Belt, resulting in occultation 
events with longer durations (jNihei et al.l 120071 ). Thus a distant occultation event would have flux 
reductions spanning severa l lightcurve points a nd requires a different analysis pipeline on the TAOS 
data from that reported in IZhang et al.l (j2008l ). Here we present our analysis of the TAOS data to 
search for objects on Sedna-like orbits. An estimate of the number of such objects detected would 
help constrain the population of bodies beyond 50 AU, thereby confronting theoretical models of 
the early Solar System, and address the puzzle of mass deflcit in the region. 
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2. Detection Algorithm 



Stellar occultations b y Solar System objects fall into the domain of Fresnel diffraction (|Born &: Wolf 



1980 



Rogues et al.l 119871 ) if the object diameter D < F, the Fresnel scale, which is given by 



where A is the observing wavelength and A is the distance from the object to the observer. For 
the TAOS system, where A ~ 600 nm, this corresponds to objects with D < 2 km at 100 AU or 
D < 7 km at 1000 AU. A high cadence system like TAOS will translate diffraction profile into 
several time-sampled lightcurve points every second as the shadow of a stellar occultation sweeps 
across a telescope. 

Figure[T] shows examples of diffraction profile for a 10 km diameter object occulting a 0.01 milli- 
arcsecond (mas) background star (V=ll AOV star) at 200 and 1000 AU, with a zero impact 
parameter (defined as the line-of-sight distance between the centers of object and background 
star). Note that we integrated every 0.2 sec with a random starting time to simulate the 5 Hz 
sampling of the TAOS observations. In order to find possible events in the TAOS lightcurves, our 
task is to detect such shallow flux reductions in consecutive lightcurve points. 

To detect occultation events by objects in Sedna-like orbits, we have applied an algorithm to 
identify consecutive, relatively shallow flux drops. The key is to recognize systematic, weak flux 
reduction with an unknown duration. For each point in the lightcurve, we defined two windows 
consisting of a set of consecutive photometric measurements centered on that point. First, a 
window {background window) of a certain size is selected, in which the median value serves as the 
local nominal fiux level. Then, at the center of the background window, a smaller window (signal 
window) is chosen that is relevant to the size and distance of the objects we want to search for 
(see Figure[2]). The sum of deviations of all data values from the local nominal value within the 
signal window then is a measure of systematic brightening or dimming in the lightcurves. A large 
cumulative deviation occurring simultaneously in multiple telescopes suggests the possibility of an 
occultation event within the signal window. The concept is similar to the equivalent width (EW) 
in spectroscopy for the strength of a spectral line relative to local continua and is similar to what 



Rogues et al.l (|2006l ) used in their work. We define the EW index as 

i+m 

'if) 



E (1) 



where (/) is the median value from all the measurements within the background window, and fj is 
the flux of the j^^ lightcurve point within the signal window {i — m,i + m} centered at i. The EW 
algorithm works because flux reductions of consecutive lightcurve points, albeit weak, accumulate 
to result in a small EW index. To search for events, these windows were centered at each point 
in the lightcurve sel^l, except the beginning and end of a lightcurve set where we cannot define a 



^We define a lightcurve set as a set of lightcurves of a single star taken synchronously with multiple telescopes. 
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complete background window, and the corresponding Wj value was calculated. 

The size of the signal window is determined by the occupation event width and the relative 
velocity between the observer and the object. The relative velocity is estimated as 



Vrel{(l),A) = V(s fcOS(/>- j (^^ " V ) ) , (2) 




where = 29.8 km/sec is the orbital speed of the Earth, Aau is the distance from the object to 
the observer in units of AU, and (j) is the opposit i on ari gle defined as angle from opposition to the 



line of sight to the background star. iNihei et al.l (|2007l ) gave an approximation of the occultation 
event width as 



H 



(2V3F(A,A))i + ' (3) 



where D is the diameter of the object, and is the angular size of the occulted star. The size of 
the signal window is given by n = sH/vrc\, in which s = 5 Hz is the TAOS sampling rate. As an 
example, u = 5 for a 0.05 mas star and 3 km diameter object at 300 AU and (j) = 30°. We have 
confirmed empirically that the choice of the width of the background window is not critical, as long 
as it is sufficiently larger than the signal window but not so large as to average over any slowly 
varying trends. A typical value is about four times the size of the signal window. 

A noticeable negative EW index indicates a significant local cumulative flux drop over the 
signal window, possibly due to an occultation event. However, with single telescope data, one 
cannot rule out the possibility that this signal has been caused by some random or other unknown 
processes. The advantage of a multiple-telescope system, like TAOS, is the ability to discriminate 
against such false positive events. 

We identify simultaneous low EW indices by ranking the EW indices in each telescope from 
the lowest to highest, that is, in a single lightcurve, the lowest EW measurement in the time series 
would have a rank r = 1, and the largest EW index would have rank r = A^p, where A'p is the total 
number of EW indices from each telescope. We then calculate the parameter 

"■=-'"n(t). w 

where r^fc is the rank of Wi from telescope k, and is the total number of telescopesH If 
the EW indices are low on each telescope at a single time point, the rank product should also be 
very low, giving rise to a large value of rj. This is illustrated in FigureO The left panel shows a 
small section of a long simulated lightcurve (27,000 lightcurve points for a typical TAOS observing 
run) with an occultation by a 10 km diameter object at 1000 AU with (f) = 30°. An occultation 



■^This is analogous to the approach used in I Zhang et al. 



20081 ). 
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was implanted at t = 200 seconds, and is not obvious at all in any of the single lightcurves. 
However, after processing with the EW algorithm on the three lightcurves, the event stands out 
clearly (Figure[3l right), demonstrating the power of multi-telescope data and our analysis pipeline. 
The values rjfc are uniformly distributed and discrete over {1, A^p ). If we make the approximation 



that the values are continuous, it can be shown (see iRicd 120071 . for example) that if there is no 
correlation between the data from different telescopes, and if the variance of the EW index in a 
filtered lightcurve remains constant throughout the duration of a data rur0, the distribution of rj 
follows the F-distribution of the form 

„iVT-l 

This is illustrated in FigureHJ which shows a histogram of ij for a single data run, along with a 
histogram of r] where one synthetic event was added to each lightcurve set and the test was repeated 
22139 times to get the distribution. While many of the synthetic events are marginal and have low 
values of ry, there is clearly a significant population of events at large values of rj. To set a cut on 
the parameter r], we note that rj depends strongly on the value of iVp, the length of the lightcurve 
set. Since, as discussed in ^ this value can vary significantly due to data runs which are shortened 
by bad weather, we can derive a parameter which would be valid for all data runs by integrating 
the tail of the distribution to calculate 



oo 

»? 

= TiFW 

In this paper we use F{r]) to select possible events. F{i]) is well-suited to this application, 
as illustrated in Figured] which shows the histogram for a data run in which there are no evident 
events, and also a histogram of F{ri) for the same data run, but with artificial events inserted into 
the same data. 

We can choose a threshold value of F that will give a few (if any) candidate events yet still 
ensure a sufficiently high event detection efficiency. We thus set a threshold of F < 10~^ , which was 
chosen empirically after analysis of a subset of data. The validity of this approach is demonstrated 
in 21 where we show that our detection algorithm gives a significant detection efficiency in the 
absence of any actual detected events. 

With further analysis we will be able, in principle, to use this approach to put a rigorous upper 
bound on the expected false positive rate. We do not do this in this paper, and the absence of 



^We define a data run as a series of tliree-telescope observations of a particular star field. 
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detected events, as shown below, legitimizes our decision. We determine our detection efficiency 
(see Section 4) by direct simulation, and thus translate our absence of detections into a formal 
upper limit on the populations of objects with orbits similar to Sedna's. 



3. Processing The TAOS Data 



A customized pipeline using apertur e photometry was d eveloped to process the TAOS image 
data. This process is described in detail in I Zhang et al.l (|2009l ). but a short summary is given below. 



At the beginning of every data run, a series of stare-m ode images is collected. These images 
are used to find the position of target stars using SE xtractor (IBertin Arnouts Ill996l ). Then these 
stars are cross-matched with USNO-Bl star catalog ([Monet et al.ll2003l ) and the position of all stars 
in the image with TAOS instrumental magnitude Mtaos < 13.5 are extracted. An aperture mask 
is then created for each star identified in the stare mode image. The aperture size for each target 
star is optimized by the signal-to-noise ration (SNR) of the lightcurves from first 1000 photometric 
measurements, and the aperture size that gives the highest SNR is chosen and used. 

The first step in analysis of the time series zipper n iode images is the subtraction of the 
sky background and the streaks from the brighter stars (see lLehner et al.ll2009l . for a description of 
zipper mode imaging and the origin of the streaks in the images). This was done by subtracting the 
median value of each column which effectively removed the background and streaks simultaneously. 
FigureO shows a subsection of a raw zipper mode images and the same image after the background 
and streak subtraction. 

The next step in the photometry process is to apply the aperture mask and add the photon 
counts. The position of the masks was adjusted at each timestamp to correct for any offset due to 
tracking errors. The offset is calculated by calculating the centroids of 20 bright stars collected at 
each epoch in the data run. 

During the photometry pipeline process, various flags were assigned to reflect known effects 
that would cause faulty photometric measurements, for example, for stars that are near or at the 
CCD image boundary or with incorrect exposure times due to system timing errors. These errors 
were all recorded and the corresponding lightcurve points were flagged and removed. 

A complete TAOS data set from 2005 February 7 to 2006 December 31, which amounts to 
~ 4.5 X 10^ three-telescope measurements was used to search for possible events by objects in 
Sedna-like orbits. Only measurements from observing runs with more than 10,000 measurements 
were included in the analysis, because short duration runs could not provide long enough baseline to 
derive significant statistics. (The duration of a normal data run is 90 minutes, however, occasional 
interruptions by system failures or bad weather would prevent the system from taking data with 
the normal duration.) Removing these data runs leaves a total of 199 data runs which amounts to 
1.7 X 10^ star-hours. A plot of the number of star-hours versus SNR for this data set is shown in 
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Figure[6l Furthermore, an a dditional 26 data ru ns were found to have significant correlations among 



the lightcurves for each star (jZhang et al.ll2008l ). Removing these runs leaves us with 173 data runs, 



corresponding to a total of 1.7 x 10^ star hours. 

Meanwhile, some of the earlier observing runs suffered from guiding problems in the beginning 
of each run. They all happened near the first 20 seconds when the guider tried to adjust the tracking 
for the first time. The common feature for this problem is a large displacement in both RA and Dec 
coordinates of every star. Hence, the first minute of data runs with these large positional offsets 
were removed. 

After all bad lightcurve points were removed, the EW filter was applied to the remaining 
lightcurves. Any rank triplets beyond our threshold of F < 10^^ are flagged as candidate events. 
Most of the candidate events identified by the EW filter turned out to be time coherent. Figured 
shows two such cases of the number of candidate events within a short period of time (ten seconds). 
For example, a run with possible high altitude clouds or cirrus would cause overwhelming flux 
reductions for most stars at the same time. Data runs exhibiting such effects are removed from 
further analysis. A total of 31 data runs were removed by this cut. On the other hand, during a run 
with low-lying moving clouds, coherent events would happen only sporadically, when the clouds 
cover some of the stars in the target field. Data that show such time dependent patterns, which 
are well characterized and understood, would not be analyzed further for event detection. 

As a final selection criterion, we required all candidate events to have at least three consecutive 
flux measurements below the nominal background level (/) used to calculate Wj. In the end, no 
candidate occultation events from objects in Sedna-like orbits were found in the first two years of 
the TAOS data. 



4. Efficiency Calculation and Event Rate 

For a given lightcurve set I, the differential occultation event rate Ri as a function of diameter 
D and distance A is 

(fRi{D,A) ^ dMD,A)vrcii^i,A)Hi{D,A) 

dDdA dDdA A2 ' ^ ^ 

where n{D,A) is the number density of objects with diameter D at distance A, Vj^ci{4>i,A) is 
the relative velocity between occulting object and observer during the data run (see Equation[2]), 
and Hi{D,A) is the event width, which is assumed to be the diameter of the first Airy ring of 
the diffraction profile as discussed in ^ (The subscript on 0/ indicates that the relative velocity 
depends on the angle between the target star and opposition during the data run, and the subscript 
on Hi indicates that the event cross section depends on the angular size of the target star.) 

With the event rate, we can now calculate the number of events expected to be seen for a given 
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star during a given data run, as 

d^iVexp,(Z),A) ^ d^Ri{D,A) 
dD dA dD dA 



EieiiD,A), (8) 



where Ei is the duration of the data run, and ei{D,A) is the efficiency for which occultations by 
objects of diameter D at distance A are detected in the hghtcurve set during the data run. 

The efficiency £i{D, A) was estimated by a simulation, in which we implanted occultation events 
in the actual light curves, and processed with the same EW filter. If we added A^^add; (-C^j A) events 
with diameter D and distance A to a lightcurve set and recover iVrec; iD,A) of them, the detection 
efficiency is simply ei{D,A) = A^rec; A)/A^add; A). Even though the impact parameter and 
epoch at which an occultation takes place affect if a particular event can be detected, by choosing 
these parameters randomly over their true distributions, we automatically average the efficiency 
over these parameters, and they can be ignored. We thus chose the impact parameters uniformly 
in the interval [—Hi/2, Hi/2], and the event epoch uniformly in the interval [to, to + Ei], where to 
is the start time of the data run. 

With £i{D,A) known, the total number of events expected to be seen in the entire data set 
can be estimated from Equation[8l That is. 



d^N,,p{D,A) _yd^Ri{D,A) 
I 



dD dA ^ dD dA 

where the sum is taken over all lightcurve sets in the TAOS data set. 



^ d^RAD.A) 

1^ wn^A ^'^hAA), (9) 



Given the extremely large number of lightcurve sets, we found we could calculate a statistically 
significant value of the detection efficiency by adding only one simulated event to each lightcurve 
set. If we add events with diameter D to a fraction wd all lightcurve sets, and events at distance 
A to a fraction of all lightcurve sets, we can write 

NsAdi{D,A) = WDWA- 

We thus choose a set of diameters D and a set of corresponding values oiwf) such that J2d = Ij 
and a second set of distances A and corresponding values of wa such that Y2a ~ ^■ 

Now there would be only one event in a given lightcurve set, so, taking into account of the 
weighting factors oi wd and wa., we have the total number of expected events, expressed in the 
fully expanded form, 



d^N,,p{D,A) _dMD,/^) 1 yV,,i{(Pi,A)Hi{D,A)Ei 

I 



^ v,,l[<Pl,A)Hi{D,/\)Ei 

2^ A2 ' 



dDdA dDdA wdwa ^ A2 

where the sum is now over all lightcurve sets where an event of diameter D and distance A is 



recovered. This expression can be simplified as (jZhang et al.ll2008l ). 



fN^D^_fniD^ 

dDdA ~ dDdA ^'') 
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where 

is the effective solid angle of the survey. 

Equation[TT] shows that the number of expected events A'exp in our survey depends on two 
factors, a model dependent factor n and a survey dependent factor Qe- This indicates that the 
number of events detected in the survey and the value of determined by the efficiency calculation 
can thus be used to place constraints on n{D, A). 



As in IZhang et alJ (|2008l ). we used the event simulator described in lNihei et alJ (120071 ) in our 
efficiency calculation. The distances, diameters, and weighting factors used are shown in Tabled! 
The weights were chosen to give preference to smaller objects where the detection efficiency is 
expected to be low in order to improve the statistical accuracy of l^c- After implanting simulated 
events, the exact same selection criteria described above were applied to find which events are 
recovered. Several example lightcurve sets where simulated events are recovered are shown in 
Figure[8l and the resulting values of Jig are shown in FigureO 

The ratio of recovered to implanted events for each value of D and A is shown in Table [2j 
Note that the algorithm did not recover any events with D < 2 km. Also note that the efficiency 
for detecting 30 km objects at 100 AU is only about 30%, indicating that a large fraction of the 
low SNR lightcurve sets used in the analysis are not useful for detecting occultation events at these 
distances. While including these stars in the analysis lowers the detection efficiency, they have no 
effect on the resulting values of Og- However, in future analysis runs we will likely set stricter cuts 
on SNR values to use our computing resources more efficiently. 

Given 0,^ and the fact that no events were found by the survey, we can set constraints on 
n{D,A). Because there are no published models of the size and distance distributions of objects 
in this region, we place constraints on n{D,A) at each value of A listed in Table [H For the size 
distribution, we assume a simple power law model of the form 



dD 'VI km, 

where q is the slope of the power law, and no is the number density of objects with D > 1 km at 
distance A. The expected number of events at distance A is thus 



A^e,p(A) = no(g-l) / ( ) ^-(^'^)^^' (1^) 

where -Dmim for each distance, is the detection limit of size below which the efficiency goes to zero. 
We picked -Dmax = 30 km because we expected the number density of larger objects to be so small 
as to provide a negligible contribution to the event rate. 

The sensitivity of the TAOS system to the size distribution is shown in FigurefTOl For q = 3 
the the survey has the maximum sensitivity for objects with D = 3 km at A = 100 AU, and the 
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Table 1. Distances, diameters, and weighting factors used in efficiency test. The weighting 
factors are normaUzed so that YId = 1 and wa = 1- 



D (km) 


Wd 


A (AU) 


Wa 


0.5 


100/785 


100 


0.2 


0.7 


100/785 


200 


0.2 


1.0 


100/785 


300 


0.2 


1.3 


100/785 


500 


0.2 


2.0 


100/785 


1000 


0.2 


3.0 


100/785 






5.0 


100/785 






10.0 


50/785 






20.0 


30/785 






30.0 


5/785 







Table 2. Distribution of synthetic events in efficiency test 



D (km) 


100 AU 


200 AU 


300 AU 


500 AU 


1000 AU 


0.5 


0/6409 


0/6325 


0/6496 


0/6361 


0/6378 


0.7 


0/6393 


0/6339 


0/6392 


0/6420 


0/6411 


1.0 


0/6339 


0/6363 


0/6466 


0/6365 


0/6453 


1.3 


0/6260 


0/6391 


0/6291 


0/6357 


0/6281 


2.0 


21/6403 


3/6432 


0/6382 


0/6285 


0/6458 


3.0 


68/6408 


23/6366 


5/6375 


1/6427 


0/6296 


5.0 


207/6302 


96/6402 


57/6357 


21/6398 


2/6318 


10.0 


291/3116 


203/3227 


172/3285 


86/3246 


29/3221 


20.0 


327/1897 


291/1927 


257/1936 


193/2010 


96/1855 


30.0 


78/295 


66/311 


69/366 


51/331 


26/310 
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diameter of maximum sensitivity decreases as q is increased. Note that the sensitivity drops off for 
larger objects, justifying our choice of -Dmax = 30 km. 

We are now ready to bring in observational results to set constraints on the Sedna-like popu- 
lation, including the number density and size distribution. The null detection in the TAOS data 
means that any model which predicts A'^cxp > 3 can be ruled out at the 95% confidence level (c.l.). 
Resulting upper limits on no as a function of q are shown in Figure[TTJ Most of our target fields 
have ecliptic latitudes \h\ < 10°, with more than 93% of the data collected in fields with \b\ < 3°, so 
the upper limits we derived on uq are valid along the ecliptic, assuming that the size distribution 
is independent of ecliptic longitude. 

We know very little about the size distribution of objects in Sedna-like orbits; Sedna is the 
only object that has been found. Given the assumed power-law size distribution, we calculated the 
upper limits on the density of objects larger than Sedna by integrating Equation ll3l 

where Dg = 1600 km is the diameter of Sedna. This upper limit on Us vs. q at 100 AU is shown in 
FiguredU Because one Sedna has been found near 100 AU, we can exclude any model that predicts 
fewer than 0.05 objects with D > Dg at the 95% c.l. That is, with the power-law size distribution, 
any size distribution which yields Ug < 1.21 x 10~^ deg~^ is excluded. An interesting consequence 
is that slopes q > 5.4 are also excluded. A larger slope means either too few large-sized objects 
(D=1600 km) to be consistent with the existence of Sedna, or too many small-sized objects (> 
1 km) to comply with the null detection by TAOS. 



Roques et al.l ()2006l ) reported the possible detection of two occultation events consistent with 
sub-km objects near 200 AU. The resulting 95% c.l. upper and lower limits on no vs. q are plotted 
in Figure[T3l along with the 95% c.l. upper limit set by TAOS survey at 200 AU. Again, if the 
assumption of a power law size distribution is valid, the interpretation of these two events as true 
occultations is consistent with the TAOS upper limit at the 95% c.l. only if q > 7.6. Note that this 
result is valid for any population of objects (e.g. Sedna-like objects, the Scattered Disk or Inner 
Oort Cloud) in this region. 



5. Conclusion 



We have analyzed the first two years (2005 February to 2006 December) of TAOS data to 
search for objects in Sedna-like orbits. With the null detection of objects in this region, we are 
able to set upper limits upon the number density of objects at different distances. Given that one 
Sedna has been found, we have also placed lower limi ts on the nurnber de nsity of objects at 100 AU. 
We also show that the candidate events reported in lRoques et al.l (|2006l ) are inconsistent with the 
TAOS data at the 95% c.l. if the size distribution at A = 200 AU has slope q < 7.6. We plan 
to modify the detection algorithm to facilitate an exact calculation of the statistical significance 
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of candidate events, and we will subsequently complete the analysis of the second TAOS data set 
(January 2007 to December 2008) and report its results in a future paper. We will also expand the 
analysis to include limits on the populations of Scattered Disk and Inner Oort Cloud objects. 
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Fig. 1. — Diffraction profiles for 10 km diameter object occulting a 0.01 mas background star 
(projected stellar diameter ~ 7 km at 1000 AU) at 200 AU and 1000 AU. Solid lines are the 

diffraction profiles, solid circles represent 0.2 sec (5 Hz) sampling light curve points, dotted lines 
indicate the duration of the 0.2 sec integration. The sampling starts at random time offset. The 
velocity of the object was calculated at an opposition angle of 30°. 
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Fig. 2. — In this simulated lightcurve, the local background level is estimated in the background 
window. After subtracting the nominal level, the summed value within the signal window, i.e. the 
EW index, is used to search for a possible occupation event. 
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Fig. 3. — The weak flux drop for consecutive lightcurve points (inset, left panel) is not readily 
discernible in the single simulated lightcurve (we plot only 3000 out of 27000 lightcurve points from 
one telescope for clarity), however, it stands out clearly as a very strong signal after processing 
simulated multi-telescope data with the EW algorithm, as shown in the right panel. 
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Fig. 4. — Thin line: distribution of rj from a data run with no occultation events. Thick hne: values 
of T] where one synthetic event was added to each lightcurve set and the test was repeated 22139 
times to get the distribution. The synthetic events are from 5 km diameter object at 300 AU. 




Fig. 5. — The Background and streak removal was done by subtracting the median value of each 
column which effectively removed the background and streaks left behind by bright stars. 
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Fig. 6. — The histogram of SNR from all lightcurves with more than 10,000 data points. 
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Fig. 7. — Histograms of the number of candidate events vs. time, with ten second binning. A 
run with possible high altitude clouds or cirrus (left) would cause frequent flux reduction to many 
stars within a short period of time. Many time coherent events could be observed throughout the 
entire run. On the other hand, in a run with moving clouds (right), only sporadic coherent events 
happened at the time when the clouds moved into the our target field. 
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Fig. 8. — During the efficiency test, synthetic events were added into TAOS hghtcurves. Examples 
of hghtcurves with the synthetic events in three TAOS hghtcurves are shown in the upper part of 

each panel and the corresponding ij values in the lower part of each panel. Some constants were 
added to second and third hghtcurves to separate them for clarity. The dotted-lines are the rj value 
for F < (a) 2 km, 200 AU, (f> = 68.8° (b) 3 km, 500 AU, = 45.9°, (c) 5 km, 1000 AU, 

(p = 34.4° (d) 10 km, 1000 AU, (p = 34.4°. 
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Fig. 9. — Effective solid angle versus size for different distances from the TAOS efficiency calcula- 
tion. 
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Fig. 10. — A plot of the the product of the power law size distribution {D /I km)^^ and ^c{D, A) 
versus diameter D at 100 AU, for several values of q. Note that D~'^ decreases with D, and 
increases with D. The product of these two terms is an indicator of the sensitivity of the TAOS 
system. 




Fig. 11. — The upper limits of number density with objects larger than 1 km for various q and A. 
Any model with no above the relevant line is ruled out at the 95% c.l. 
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Fig. 12. — Upper limit on as a function of q at 100 AU. The diagonal line is the 95% upper 
limit set by the TAOS survey at 100 AU. The right axis shows the corresponding upper limit on 
the total expected number of objects larger than Sedna (A^s) in the whole sky,assuming an isotropic 
distribution. Given that one Sedna actually exists, we can also set a lower limit on the surface 
density. If our assumption of a power-law model is correct, we can exclude size distributions where 
As < 0.05 (dotted line) at the 95% c.l., which corresponds to a value of rig < 1.2 x 10~^ deg~^. 
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Fig. 13. — The 9 5% upper and lower l imits on number density were calculated based on two ~ 600 m 
objects found bv lRoques et al.l (j2006l ) and q from 3.0 to 8.0, along with the 95% upper limit set by 
TAOS survey at 200 AU. 



